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Abstract 

r) , This article addresses two topics of significant mathematical and practical interest in the 

theory of kernel approximation: the existence of local and stable bases and the Lp-boundedness 
of the least squares operator. The latter is an analogue of the classical problem in univariate 
2 ' spline theory, known there as the "de Boor conjecture" . A corollary of this work is that for 

appropriate kernels the least squares projector provides universal near-best approximations for 
functions / € Lp, 1 < p < oo. 

(N' 

> 

(N '. 1 Introduction 

(N 

'^ I The purpose of this article is two-fold: to provide operator norm estimates for least squares ap- 

^-^ ■ proximation using linear combinations of translates of kernels, and to demonstrate that there are 



^ 



o 



c^ 



local, stable bases for the basic spaces of kernel approximation. We consider the problem of finding 
the best L2 approximant from the finite dimensional space 



r> I 5(k, H) := spanK(-,^), 



where H C M, M is a d— dimensional, compact, complete, smooth Riemannian manifold, and 
At: M ^ M is a positive definite kernel (we make more assumptions about it, anon). In pursuing 
this result, we identify a basis (w^)ges for 5'(k, H) having uniformly bounded condition numbers. 
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A large body of the literature has been devoted to the interpolation problem: given /, find /=/ € 
S{k, E) by enforcing the condition 

iI^f)\E = fk- (1.1) 

The cornerstone of the theory is that the linear operator 1= is the orthogonal projection on a 
certain reproducing kernel Hilbert space [Nf^ (called the native space) to its subspace S{k,E). In 
other words, /=/ is the best approximant to / from S{k, H) in X^. 

Recently, the authors of [8j have proved that, under basic geometric conditions on the sets H (a 
condition defined in the following section known as quasiuniformity, which we assume throughout 
this article) the L^o norm of the interpolation operator /=/, known as the Lebesgue constant L, 
is bounded independent of H. I.e., the interpolation operators I-,, originally bounded by 1 in 
the !Nk — )• Kk norm, are projectors from C(M) to S{k, E) with Loo — > -^oo norms bounded by a 
constant independent of the dimension of S(«;, H). An immediate corollary is that for / G C(M), 
interpolation is near best approximation, as is shown in the following standard inequality: 

11/ - hfWoo < (1 + £)distoo(/, S{k, E)). 

To prove this result, a local basis for S{k, E) was utilized. In the first part of this article, we 
investigate the Lp stability of this basis: i.e., we give conditions on the kernel sufficient to guarantee 
norm equivalence of the Lp norm of a function in S{k, E) and the sequence norm of its coefficients 
(relative to the local basis). 

Our second goal is to study the least squares projector Th from from Lp(M), 1 < p < oo to S{k, E). 
For each / G Lp(M), we define Ts/ to be the unique element in S{k, E) such that 

/ (fix) - {T^f) {x))g{x)dx = 0, 5 G S{k, E). 

Jm 

When p = 2, T^f is the best L2 approximation to / G L2 from S{k, E). Because it is a projector, 

\\f-Ts{f)\\p<{l + \\mp)distp{f,S{K,E)). 

We show that having a local and stable bases is sufficient to bound uniformly the norm of the least 
square projector (independent of #H, similar to the case of the Lebesgue constant). In particular, 
this guarantees that least squares projection is near best approximation in Lp, 1 < p < 00. 

The development here is reminiscent of an episode in the history of the univariate B-spline theory, 
dubbed the "de Boor conjecture" [l], which stated that the L^q— norm of the least square projector 
onto the space of splines of order k can be bounded independently of the knot sequence. Douglas, 
Dupont, and Wahlbin |lj established the Loo boundedness of the least square projector onto the 
space of splines of order k {k > 1) with knot sequence t := (tj) in terms of the mesh ratio: 

Mt := max(tj4.i - tj)/min(tj+i - tj). 

i i 

A more stable upper bound for the projector in terms of the smaller global mesh ratio 

Ml := max(tj+A; - tj)/min(tj+fc - U). 

i i 

was given by de Boor \2\- Much later, Shadrin [19] gave an affirmative answer to the conjecture. 



Although the manifold setting is rather abstract, conducting our investigation in this generality is 
sufficient to cover many practical applications, including approximation on spheres (S*^), projective 
spaces (P'^), and many matrix groups. On the other hand, little is required in the way of geometric 
sophistication; most statements hold for compact metric spaces, and the focus on manifolds comes 
mainly because the examples of local kernel bases from [5] have been developed in this setting. 

We organize this article as follows. In Section [2l we make necessary preparation for proving the 
main results. In Section [3l we establish upper and lower condition numbers for the Lagrange basis 
under some simple assumptions. Section [J] investigates the problem of continuous least squares 
minimization. We use the condition numbers obtain in Section [3] to estimate the spectrum of the 
Gram matrix for this problem and bound the ioo — >■ ^oo norm of its inverse: we show that this 
matrix is boundedly invertible when the basis is local and stable. In Section [5l we state the main 
results for Lp stability of the L2 projector, the L2 projector is Lp stable when there is a local, 
stable basis. We then observe that this holds for the spaces generated by the kernels introduced 
in [8]. An example is given for the 2-sphere where precise approximation orders are known. A few 
conjectures are discussed in Section O while the Appendix houses some technical proofs about the 
kernels introduced in [8]. 



2 Background and Notation 

2.1 Positive definite kernels 

The motivation for this article was to consider least squares approximation from S{k, H), with k a 
positive definite kernel and H a finite subset of M. In this case, a basis can be obtained directly 
from the kernels, v^ = k(-,^), (^ € H and, until recently, it may have seemed natural to consider 
this basis for most computational problems (the reproducing property of k: 

{f{-),K{;x))=f{x), /EX«, XGM, 

perhaps motivated this natural instinct). However, such bases are neither local nor stable, which 
results in analytical and numerical difficulties. In particular, the global nature of the kernel leads 
to a full Gram matrix, G = {{!<'{■,£,), k{-,C))) c i-(z-^ ^^"^ little can be said about the norm of its 
inverse. 

Even in the case of the compactly supported Wendland functions, as the number of centers increases 
the bands of associated collocation or Gram matrices become ever larger. Instead, we wish to use 
a basis that is local, and that scales according to the density of the centers. 

Recently, in [8j, this has been demonstrated for Lagrange functions associated with interpolation 
using certain kernels. For a general compact, Riemannian d— dimensional manifold M, a family of 
kernels Km,M, for integers m > d/2 were considered. These kernels have as their native space ^Nf^^^ 
the Sobolev space W^{M). 



2.2 Manifolds 

Throughout this paper, M denotes a compact, d-dimensional Riemannian manifold. We choose to 
carry out our investigation in the manifold setting in order to make the results applicable to a broad 
spectrum of practical problems, including approximation on spheres (S*^). However, familiarity with 
Riemannian manifolds is not necessary for a thorough understanding of the arguments we present 
here. For the definitions, theorems and proofs in this article, the salient facts are that M is a 
compact metric space, and a few more conditions we now make clear. 

We denote the ball in M centered at x having radius r by B{x,r). Given a finite set H C M, we 
define its mesh norm (or fill distance) h and the separation radius q to be: 

/i := sup dist(x, H) and q '■= - inf dist(.^,C)- (2-1) 

The mesh norm measures the density of H in M, the separation radius determines the spacing of 
S. The mesh ratio p := h/q measures the uniformity of the distribution of H in M. We say that 
the point set H is quasi-uniformly distributed, or simply that H is quasi-uniform if H belongs to a 
class of finite subsets with mesh ratio bounded by a constant poO 

The manifold is endowed with a measure, and we indicate the measure of subsets fi C M by vol(ri). 
The integral, and the Lp spaces for 1 < p < cx), are defined with respect to this measure. The 
embeddings 

C(M) C Lp(M) for 1 < p < oo and Lp(M) C Lg(M) for 1 < g < p < oo 

hold. In addition, L2 is a Hilbert space equipped with the inner product (•,•): {f,g)^ {f->9) 

The final properties of the manifold, in force throughout the article, are the following. 

Manifold Properties. There exist constants < aMj-f^Mj'^M < 00 and < rM < co for which 
the following hold: 

1. for all X &M. and all r < rM , 

"'^ <vol{B{x,r)) <ujMr'^; 



2. /or S C M, 2; e M and R>q, 

#{En B{x,R)) <KM{R/qf. 

We note that the preceding holds when rM is the injectivity radius, which bounds the radii of 
balls in which local geometric properties of the manifold are similar to those of M'^o Its technical 
definition is unimportant for our purposes here, although it appears in the statement of a key 
condition (Assumption 13.4]) . as well as in the definition of the previous constants. 



^Naturally, every finite point set has a bounded mesh ratio. Our results hold for point sets from a (not explicitly 
defined). Inclusion in this family means that the mesh ratio is bounded by an unspecified constant whose size 
influences various constants in our results. 

^This is the radius of the largest ball in R'' for which the exponential map (centered at any point of the manifold) 
is a diffeomorphism between the ball and the corresponding neighborhood in the manifold. 



3 Condition of the Lagrange Basis 

In this section we show that a Lagrange basis is stable if it is local and satisfies an equicontinuity 
condition. 

Let {v^)^^:= be a linearly independent family of continuous functions on M, indexed by a finite set 
H C M. In this section we present conditions sufficient to establish an equivalence between the 
Lp norm of s € X := spamg^f^ and the ip norm of its normalized coefficients. Thus, we seek 
constants < ci < C2, called condition numbers, so that for s = X^fes ^C^?' ^^^ t^° sided estimate 

ci^^/^Ei^^n'^'^ iiE^c-diL.(M) <c2<z'^/nEi^d^)'^' (3.1) 

holds true, with the usual modification, max^g= |^^| ~ || X^fg= ^^f^||oo, when p = oo. In particular 
the basis to be considered will be the Lagrange basis (see Assumption 13. 3p . One observes that, 
as expected, the condition numbers depend in a controlled way on p, as is the case with B-spline 
bases and wavelet bases, and this dependence on p (encoded in the factors q'^'P) can be suppressed 
by taking an appropriate renormalization of the basis (see Remark 13.21 below) . 

Our ultimate goal is to prove bounds (j3.ip for certain kernel spaces X = S{n, H) with condition 
numbers that are independent of t^H. However, the development of this section does not make 
use of the fact that X is a kernel space (and there may be instances where one would not wish to 
assume this). Instead, we prefer to state our results for a general space, X, in terms of a few basic 
assumptions on the the basis functions v^, and then to show that such assumptions hold for the 
kernel spaces we consider. 



Remark 3.1. An immediate consequence of 113. 1\) is the Nikolskii inequality, which states that for 
1 < p < q < oo, 

||s||q < — q ^p'q'WsWp for all s ^ X. (3-2) 

ci 

In case X = S{k,E), the penalty q ^p i'^ is in line with other examples of Nikolskii 's inequality, 
since dim(S'(K,H)) ~ q^'^ when k is positive definite and H is quasiuniform. 

Remark 3.2. We often write s & X as a renormalized series s(x) = X^fgs ^p,?^p,5(^)' where 
Vp^^{x) := q-'^/Py^ix) and Ap^^ := q'^lvA^. 

In establishing condition numbers, we employ three basic assumptions on the family of functions 
(f^)gg=. We assume the basis consists of the Lagrange functions, the family is fairly localized 
(exhibiting rapid decay away from a center), and the basis functions are Holder continuous, with a 
uniform seminorm that scales appropriately with the centers S. We now formally state these. 

Assumption 3.3 (Lagrange Basis). The basis (u^)ggH is the Lagrange basis. I.e., 

A consequence of this assumption is that coefficients are simply the sampled function values. Thus, 
the inequality we are after can be recast as a sampling inequality 

Ci||s|h||^p(H) < q''^^^ \\4lp{M) < C2||s|h||^^(H)- (3.3) 



Literature abounds in the establishment of similar inequalities. In particular, if the functions 
considered are polynomials, then these are often referred to as Marcinkiewicz-Zygmund inequalities 
([llj. [5J, |18j); if the functions considered are entire functions of exponential type, then these are 
referred to as Plancherel-Polya type inequalities ([E], [H], [E]). 

The second assumption is a localization assumption, and is stated for a general family of functions 
(not necessarily a basis) indexed by the centers .^ € S. 

Assumption 3.4 (Localization). For a family of functions {v^)^^^, there exist positive constants 
V and Ci for which each individual function v^: M — )• C satisfies 

I .^1/^ / min(dist(x,0,rM)\ 
|u^(x)| < Giexp —u . 



Under the assumption that the set H is quasi-uniform, the following is equivalent to Assumption 

I .^|,^ I / min(dist(x,g),rM) \ 

\vi^{x)\ < Ci exp I -u '- I . (3.4) 

with v'p = V. We will use them interchangeably according to convenience. 

Similar bounds have been established in [8] for the Lagrange functions obtained from a large class 
of kernels. These are the kernels Km,M) "i > d/2, the reproducing kernels for the Sobolev space 
Ty^(M) endowed with a certain inner-producto An immediate corollary is that the Lebesgue 
constant sup^g^^Ag^ 1x^(3^)1 is bounded in terms of Ci, v, p and M (cf. [8l Theorem 4.6]). 

The final assumption is one of equicontinuity. 

Assumption 3.5 (Holder Continuity). There is < e < 1 so that for any quasiuniform set H C M 
with separation distance < q, 



\vi:{x) -wg(2/)| < C2 



dist(x,y) 



q 

with constant C2 depending only on e, M and the mesh ratio p = h/q (but not on ^^HJ. 

This final assumption is closely related to Bernstein inequalities, which have been investigated 
recently on spheres in [12], and have been demonstrated for a very large class of kernels. Namely, 
for all zonal kernels that are a perturbation of the fundamental solution of an elliptic partial 
differential operator by a convolution of the solution with an Li function (this is |121 Theorem 
5.5]). We remark that [SI Section 5] demonstrates that the kernel K2S2 of order 2 on S^ is such a 
perturbation. Using entirely different techniques, we will demonstrate that all kernels of the form 
Km M satisfy Assumption 13.51 in Lemma 17.31 in the appendix. 



^More precisely, they are the reproducing kernels for the Sobolev spaces W^™(M) endowed with the inner product 
(u, v)m = SfcLo /m^^'°'^' V°i')a;d2;, where V is the covariant derivative on M, which maps tensor fields of rank j to 
tensor fields of rank j ' + 1. Although this is of fundamental importance to establishing the localness of the Lagrange 
basis (and for bounding the Lebesgue constant) it is a technicality which does not directly play a major role here. 
We direct the interested reader to [8] Section 2] for details. 



3.1 Upper Bounds 

We depart temporarily from using Lagrange functions to prove the upper bound, since this prop- 
erty is not necessary. We treat famihes of functions that merely satisfy the localization property, 
Assumption 13.41 with centers that are quasiuniform. 

Proposition 3.6. For quasiuniform centers H C M and a family of functions (f^)5eH satisfying 
Assumption \3.4\ we write s = X^f es ^C^? ~ SfeH ^p,?^p,€- There exists a constant C2, depending 
only on M, p and the constants appearing in Assumption \3.4\ so that, 



\s\\p < C2IIA 



kp(S)- 



Proof. For p = 00, the upper bound follows by a standard argument, decomposing the sum en 
annuli. This is the argument used to bound the Lebesgue constant in the proof of [H Theorem 4.6] 
(in fact, if v^ is the Lagrange basis, the upper bound we seek is exactly the Lebesgue constant). 
Thus we have L for which 

||s|1l^(M) < £||s|s||^^(H) = ^Poo,-||i'^(E), 

For p = 1, we have 

Is(x)|dx < C||Ai^.||^j(='). 

Here we have used the fact that Hfi^^Hi < C for some constant C depending only on the manifold 
M and the constants in Assumption 13.41 A standard application of operator interpolation proves 
the other cases. The constant C2 can be taken to be C^'P£^~^'^. D 

Note that in the proof above we did not use Assumption 13.41 in its full force. In fact, the following 
much less stringent decay condition for the family of functions 

Kc(2;)| <g"''/^(i + dist(x,0/g)~^ /i>d, 

is sufficient for Proposition 13.61 to hold true. 



3.2 Lower bounds 

The lower bounds in this case are sampling inequalities, and the lower bound in the L^q case follows 
immediately with constant 1. The main challenge is to establish the Lp lower comparison, and, in 
contrast with the upper bound, there is no convenient way to make use of Riesz-Thorin or other 
interpolation results. In our proof, it is not just the decay of the basis elements, but also their 
nature as Lagrange functions and their equicontinuity that is important. 

Proposition 3.7. Under Assumptions \3. 31 \3.4\ and \3.5l there exist constants ci > and qo > 0, 

so that for q < qo, for 1 < p < 00 and for all s G X, 

ci\\s\sh,(s) < q^'^^'^Mp. 

The constants ci and qo depend only on p, the constants appearing in Assumptions \3.4\ and |g.5l 
(namely Ci, v, C2, e), and on M (specifically on the constants appearing in Assumption[27i 



Proof. The L^o case follows imniediately with constant 1. We will therefore assume in the rest of 
the proof that 1 < p < oo. 

We estimate the Lp norm of the function by integrals over a union of disjoint balls {-B(C, 7^) : ^ G H} 
with 7 < 1: 






dx. 



The Holder continuity assumption implies that, for 7 > sufficiently small |wg(a;)| > ^ in a 
neighborhood B{^,^q) (the choice of 7*^ < 1/(3C2) does the trick). Consequently, we have, for each 
fixed ^ € H, that 

/ \v^{x)\Pdx> min \v^{x)\Px {vol{B{^,jq)))> (I) auiiqf. 

We now use a diagonal dominance argument: we show that, for each fixed ^ G H, the integral 
/-BCf -Ya) I ^cge. ^iO'^ci^)\'''dx is dominated by the "principal term" Jb(c^„\ \s{^)v^{x)\Pdx. To show 
this, we estimate the contribution of the off-center terms. In short, we have the estimate 



E 

«6S 



/ \Y,s{C)vdx)'dx>Y,l2''^ [ 



\s{Ov^{x)\Pdx - / \Y.s{C)v^{x) 



dx 



This is a consequence of the following inequality which we use at times throughout the paper, and 
state here formally: 



E' 



<-^-^E 



(3.5) 



j=i 



In this case we have applied (13. Sp with n = 2, oi = Ylcer. k(C)f^c(^)l ^"^^ 02 = — YliC&'S\£ l^(C)^c(^)l' 
which after a manipulation gives |ai|^ > 2^~'P{ai + 02)^ — |a2|^- Combining this with the above 
estimate of the principal term gives 



\s\\i>Y. 



3P 



{i<iTM 



y]s(c)vc(3 



dx 



(3.6) 



By shrinking 7, the sum Ylf&'S.lBieig) I Z^c^^f *(C)^c(^)l^*^^ ^^^ '-'^ made as small as desired. In 
particular as shown in Lemma 13.81 below, the off-center sum can be made less than one half of 
the principal part of (j3.6p . X]fg= ~^{l(lY\^i,\^ ^ ^'^d the choice of 7 depends only on M and the 
constants from Assumptions 13.41 and 13.51 and not all on p. As a consequence 






Thus 



ci = min {aml'^Y'^/^ 
pe[i,oo] 

= min(aM7'^/3, !)• 



D 



Lemma 3.8. There exist two constants qo and R, depending only on the manifold 
constants from A ssumptions \ 3. ^ and \3.l^ such that the inequality 



and the 






holds for all s ^ X and all p G [1, oo) when q < qo and 'j < R. 
Proof. Given 7<l,r>l,g< F/pm and ^ S H, we write 

/ I E ^(c)^c(^) ''dx < m + 11^ + iii^), 

where, using (j3.5p with n = 3, we have 



/, := 3^ 



-1 



-1 



Ih := 3P 



Ilh := 3^-^ 



B{(rfQ) 



B{^m) 



B{S.aq) 



E <Ovd 



dist(C,C)<rg 



dx 



E dOvd^ 

rg<dist(C,C)<rM 

y s{C)v(^{x) dx. 

dist(C,?)>rM 



dx 



We will bound the three sums ^f^^I^, ^f^-I^, and ^c^^IH^ so that each is less than 
To obtain the desired estimate, we adopt the following strategy: 



1. we first choose F sufficiently large to bound Yl^h — ^'(F ^) x (7^)'^X]feH l'5(OI^' 

2. we then choose q sufficiently small to bound ^ ///^ = o(g) x (7^')'^ X^fes k(OI^) at the same 
time keeping q < rRi/F, 

3. we finally show that X^/g < {"^^ i^YilQ) z2$e'E\^iO\^ where 3<C is a constant depending 
only on M and the constants from Assumptions 13.41 and 13.5^ which allows us to choose 7. 

The order in which we carry this out is important, because both q and 7 depend on F, and we wish 
to be able to choose them independently of p. 

Step 1, estimating II: The quantity 11^ is subdivided into (comparable to) | logg|-many terms, 
each term being a sum over a dyadic region: 



Ofc := ftkiO := {C e H I T2''q < dist(C, C) < Tl'+'q}, k = 0,1, . . . ,Nq, 



in which 2^? ~ j^. This means that, for Mk := JBiiaq) ^C&^k ^^O^d 



dx, 



N„ 



fe=0 -^^(S,!?) 



^ s{C)vc_{x) 



A^o 



dx = (3)P-i ^ 2(P-i)('=+i)Mfc, 



A:=0 



where the above inequahty follows from Yll=i^j — Y^i=i'^^ W3\^ ■ This follows easily from 
repeated application of (13.51) with n = 2. We note that (I3.5p could be applied directly to estimate 
//^, but at a cost of incurring a factor Nq~ , which depends on q. 

We now estimate the contribution from each M^ , the portion of 11^ coming from the dyadic interval 
Qk- By using the inequality | Yll=i ^j\^ — n^~^ Yl Wj\^^ '^^ have 



< idtVikY'^ X max II va 11^ 

< (i^M(2'=+'r)'^)''"' X (Cf exp(-i/pr2'^ 



CGSlfc " "'-f'°o(-S(5.7'?)) 



X vo. 



CeOfc 



^miigf) X j; |KC)r. (3.7) 



In the final line, we have used the estimate #f^fc < -fCM(2 T) , and Assumption 13.41 namely that 
for C G ^fc and X G B{^,jq), we have |f;^(x)| < Ci exp(-iyr2'=). Multiplying by 3P-i2(p-i)(fc+i) and 
summing from to A^^, we obtain (after rearranging some terms) 



Nr, 



Ik < CiUMmMCi2''+y-\jqf ^(2'=('^+i)r'^)P-iexp(-z.pr2* 



, fc=0 



E i«(or 



(3.8) 



We introduce the notation JCjj := 3KmCi2'^'^^ . Summing the 11^ over ^ € H and using the 
inequality #{S,- C G ^fc(C)} ^ A'M(2'^+"'^r)'^ (this follows from the fact that, for a given ( & E, 
#{^: C G ^fc('?)} = #^fe(C)i since the condition defining these sets is symmetric in ^ and () we 
obtain, 



A^o 



E^^e < C,u;M%'u\iqr Y.(^'^'~''^^'r' eM-i^pr2')Y^ 



€es 



, fc=o 



?GS 



E i«(c)r 



CGHfc 



< CiojMX'ijdiq)'' E(2'^''"''^r'^r"'exp(-j^pr2'=)(A:M(2'=+^r)'^) 



, fc=0 



Ei«(or 



C6S 



< ^^p (79) 



d / ^^ 



"// 



rp 



^^2^P)(d+i)Pg^p(_^^P2fc) 



.fc=0 



Ei^(c)r 

C6H 



In the final inequahty, we have used the fact that 2'^''+^^'^ < 2'=('^+i) x 2'^+^ and that T'^p = ^^^^f^. 
After some minor cleaning up, we have 



^Jk < "^OC^ri (r'^eM-^pr) + 1^ l^xpi-upry^+'^P-'dr^ (jq) 



du IIP 
^lsll£p(H)- 



10 



This estimate holds because r i— )• rP("+i) exp(— z^pr) is decreasing for r > {d + l)/z^. Because we 
would like to bound this final expression by ^fr(7'7)'^||s|=||£ p-,) it suffices to take 

max('r*exp(-i.pr),^ H exp{-i^pr)r^'^+^^P-^dr^ < ^{9CiKm2'^+Y''- 

As we demonstrate in Lemma [3.91 using e~^ = -^9CiKm'2'^~^^ , this can be accomplished for a 
choice of T depending only on M and the constants from Assumption 

Step 2, estimating ///: Mimicking the estimate (j3.7p we have 



iii^ < 3P-1 (#{c G s I dist(C,6 > m})"-' X E / \siOvdx)rdx 



<3P'^{#Ef-^x max ||t.c||^ ,„,, ,, 



X vo 



dist{C,0>rM 



It is easy to estimate the cardinality of H, since vol(M) > X^cg^ ■^ol(i?(^, (7)). By using Assumption 
331 to estimate max|^idist(c,$)>rM} ll^cllL^(B(5,7g))' "^^ obtain 



III. < 3^ M -r- Cf exp -up — ojMiiq) 

= Cia;MX?7/g-^(^-i)exp ("^^Py) (t?)' 
with JCiii := 3Civol(M)/aM- Summing over .^ G H, we have 

Y^IIh < Cia;MX?7/g-^(P-i) exp (-up'-^) (79)'^ ^ 



E i«(oi^ 



dist(C,C)>rM 



dist(C,0>i'M 



E i^(OP 



dist{C,C)>rM 



< C,c.m3C?7/.-^(-^) exp f-.p^) (7,)'^ E f ^) 1^(01 
= -Y^iiiQ ^ exp I -i/|j— 1 (7^) 



E 1^(01^ 



We wish to ensure that J2fes^^h — 5§^(')'^)'^X]feH l"*(Ol^i so it suffices to take q < qo, where 
Qo ^ fMl/r is the largest value for which the condition < q < qo implies q~'^ exp {—urM/q) < 
^(9C7i vol(M))-i. 



Step 3, estimating /: To estimate /^, we make use of bounds on V(^ derived from the Holder 
continuity assumption and the fact that the Lagrange function vanishes at S,, when C 7^ C; to get 

k<^'"' [ (i^Mr'^r'(C27T E iKordx 



dist(c,C)<rg 



dist(C,C)<rq 
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In the first inequality, we use the bound 

l^'c^la;)! < C2Y ■: for x satisfying dist(x,.^) < 7(7. 
and the estimate on the number of centers 

#(Hni3(e,rg)) <KmT^. 

The third inequality follows from the simple estimate \o\{B{^,^q)) < ujMijQ) ■ 
Summing over ^ G H, we obtain: 

where % := 3KmC2- The final estimate results by exchanging the two summations, and employing 
the fact that #{C G H: dist(C,0 < ^q} < -^Mr"*- 

Thus, in order to force E^eB^S ^ l^i^lf\\''^\sfe (=)' ^^^ ^^°^^^ o^ T' < ^(9i^MC2r'^)-i suffices. 

This completes the proof of the lemma. D 

Lemma 3.9. Let i^ > 0, d > 1 . For every e > there is Tq > {d + l)/i^, (depending on u, d, e) so 
that for r > To and for all p € [1, 00), 

max (r'^Pexp(-z/pr), -^ / exp(-z/pr)r('^+i)P-Mr J < e^. 

Proof. It suffices to establish that 

2 ''°° 



— I exp(-i/pr)r('^+i)P-i dr < e^ 



2(d+l) - 



for an appropriate Fq and F > Fq. To see this, assume Fq > max(e, -^^ -) which implies 

r < — e**' " whenever r > Fq. 
e 

Hence 

— / eyip{-ijpr)r^'^^^^P^^ dx 



o /"OO 

< — / exp(-i/pr) 



-| (d+i)p-i 



-e 



dr 



dp 



< 't) tS^I -^■'- 



± r . ^^ . e-P'^r/2_ 



e y TopueP 
The desired conclusion now follows easily. D 
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3.3 Stability of the Lagrange Basis 

From Propositions 13.6 1 and 13.7] we now have: 

Theorem 3.10. Under Assumptions {3^ \3-4\ and lS.Sl there exist constants < ci < C2, so that 

Cl\\Ap,-\\ip(S) < I|s||Lp(M) < C2\\Ap^.\\ip{S)- 

holds for all s = J2fes ^C^C ^ ■^■ 

It follows that the Lagrange basis considered in [8] is stable. 

Corollary 3.11. Let Wl be a d- dimensional, compact, Riemannian manifold, let m > d/2 and let 

Km,M be one of the kernels considered in f^. If 3 C M has mesh ratio p, then the associated Lagrange 
basis {xs,)(eE, Xs, = X^ oc'^™.m(') C); satisfies /i3. 1\) for all s € S{Km,M), with ci,C2 depending only 
on p, m and M . 

Proof. The Lagrange functions satisfy Assumption 13.41 by ||8j Proposition 4.5]. It remains to show 
that they are uniformly Holder continuous. This is shown in Lemma 17.31 which we postpone until 
the appendix. D 

4 The Loo-norm of the L2-projector 

We now discuss the problem of best L2 approximation from finite subspaces of M, an elementary 
problem, but one which we consider for the sake of clarity. Given a family (indexed over the 
set H, with #H < 00) of continuous functions (f^)gg= on M, we denote the 'synthesis map' by 
V: C" -^ X := span£g5i;g C C(M). This is an operator which maps sequences of coefficients to 
functions in X via 

The corresponding 'analysis map' 



V*: Li(M) ^C-: s^ {{s,v^)) 



C6H' 



is the formal adjoint of the synthesis map. Since each v^ is in Li(M) n Loo(M), it follows that V* 
is a bounded operator on each Lp space. 

If V is 1-1, it follows that the L2 projector on X is the operator Ts = V{y*V)~^V* , and it, too, 
is bounded on each Lp(M),l < p < 00. One can write T'^if) = X^tgH"^?^?' f°^ / ^ -^P' with 
coefficients c^ (uniquely) determined by the Gram matrix G = V*V and inner products with /, 

G(ce)cGH = ^*/ ^ {T^{f),v^) = {f,v^), eeX 
A consequence of this fact is that (T^/ — f) J- X and 

||/-T=/||2<||/-S||2 
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for all s ^ X. 

We make special note of the fact that the L2 projector is independent of the choice of the basis. As 
such, the formal dependence of T=; on H should seem out of place. We adopt this notation because 
when X = S{k, H), the index set S C M is contained in M, and the spaces (not just the bases) are 
naturally indexed by H. 

The norm of T5 on Lqo matches the operator norm on Li, and this quantity controls the Lp operator 
norms for all 1 < p < 00. Indeed, 

||7s|loo= sup sup \{Tsf,g)\= sup sup \{Tsg, f)\ = \\Ts\\i. 

ll/l|oo=l ||5||l = l ||/||oo=l ||9l|l=l 

We will establish an upper bound for ||T=||oo by separately estimating the following three norms: 



\V\\e _,r , \\{V*vy^\\e ^p , IIV"*! 



Lr^^l. 



00 ^■'-oo 



The challenging part lies in the process of estimating ||(^*1^) ^Ikoo^^oo- Under basic assumptions 
on the basis functions, we can bound the product of the remaining two operator norms ||V^||oo x 



|^*||oo, with little effort. 



4.1 The £00— boundedness of the inverse of the Gramian 

Using the basis (^2,g)^GH) renormalized as in Remark 13.21 the L2 projector T-= has Gram matrix 
G := ((i^2,$)'y2,c))f,C- '^^^ equation 

Eg,ce5QgG(g,C)ac _ II Eggs^g^a.^llj 
EceHl«Cp ~ l|a|l?2(s) 

shows that the Riesz bounds associated with this basis determine the numerical range of the 
Gramian. By estimate (|3.ip . the eigenvalues of G are bounded below and above by cf and c^. 

We now investigate G and its inverse. For a local basis, as in Assumption 13.41 it is possible to 
estimate the size of an entry by its distance from the diagonal. For ^, C G H sufficiently far apart, 
we write Jli := {j; G M | dist(x,^) < dist(x,C)} and il2 := {x G M | dist(x, C) < dist(x,0}- Noting 
that both dist(C,f^i) and dist(^,J72) are no less than dist(^,C)/2, (I3.4|l implies 



|G(e,c)i 



V2,£,ix)v2^(ix)dx 



< Cg exp 



-II- 



< / \voo,i{x)\\vi^(;{x)\dx + / \voo.(:{x)\\vi^^{x)\dx 
min(dist(^, C),rM 



with // := 1^/2. Here Cg is a constant depending only on the manifold M and the constants in 
Assumption 13. 4[ 

To show the ^00— boundedness of the inverse of the Gramian, we adapt an argument developed 
by Demko, Moss and Smith [3]. Their results primarily apply to bandlimited matriceqj, and our 

^ Demko, Moss and Smith also addressed a few exceptional cases in which the matrices involved are full. However, 
the discussion there does not seem to imply the situation we are dealing with here. 
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Gramians are not such matrices. To use their argument, we decompose the Gramian as a sum of 
a bandhmited matrix B and a neghgible matrix R. To be precise, we write 

G = B + R = B(ld + B^R), 

in which B(^, C) = for dist(C, > ^Q and R(C, C) = for dist(C, < Tg, where T > 1 is a cutoff 
parameter we need to administer with care in our proof. The main challenge in the process of 
bounding ||G"^||oo = ||(Id + B~^R)~^B~^||oo is to show that ||B"iR||oo < 1. The key idea is that 
we bound the norm of R by the integral J^ r e'^'^dr, which decays exponentially with T whereas 
the norm of B~^ is comparable to V^. 

Proposition 4.1. Suppose the basis ('f^g)^gs is stable and local. I.e., it satisfies 113. 1\) and Assump- 
tion's^ Then there is a constant R > so that for all H with q < R, the inverse of Gram matrix 
G = ((t'2,5)f^2, ())§,( has Lqo norm bounded by a constant depending only on M and the constants 
from /i3.1\) and Assumptions \3.4\ In fact, we have 



-111 

lloo 



< 2er'^ 



with constants C and T (depending only on M, Ci, v and c\ and C2) as defined in ^4-^ o,nd J^. j[ j 
below. 

Proof. By picking F sufficiently large and q sufficiently small, ||R||oo < sup^ X^distff c)>rg |G('^iC)l 
can be made as small as desired. Indeed, we have 

00 

iiRiioo<sup Y. iG(e,c)i<supj; Y. iG(e,c)i 

We have defined Hk{C,) := {^ € H : kVq < d{^, () < (A;+ l)rg} to be the set of centers in an annulus 
around C of radius kTq. For small A;,^ G -f^fc(C)) ensures that |G(^,C)| < CGj2feH (c)^^P(~l^^^)^ 
while for large k, (say k > k*) ^ G -f^fc(C) implies |G(^,C)| < CGY^^(,Hk{0^'^P(~l^^^/l)- ^^ 
symmetry, the same inequality holds for ||R||i. Since the cardinality of each set -^^^(C) is less than 

Kmik + 1)'^T'^, while the cardinality of UA;>fc* ^k{C) C H is less than ( JEE^ — 1 J ^ t^q have 
max(||R||i,||R|U) <Cg/^m | F-^ f^ik + l)-^ exp(-/xfcF) + (" ^^^^^W \ expi-firu/q)] 

< X (F'^exp(-^F) + i /""r^e-^Mr + r^^^^!^^') exp{-fim/q)] • (4.1) 

By choosing F sufficiently large and then q sufficiently small, we can impose the condition that the 
last expression is less than 

imin(c?,e-iF-'^), 
where S is defined in Equation (j4.3p . 

It follows that ||R||p is less than gcf for every 1 < p < 00. In particular, this holds when p = 2, 
which implies that cr(B), the spectrum of B, is contained in the interval [^cf, |c|] C [^cf,2c2]. 
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We are now ready to use Demko, Moss and Smith's argument. It is possible to approximate the 
univariate function - on [2C1, 2C2] using pn, an algebraic polynomial of degree n with uniform error 



max 



■^^19^1' 2 



Pn{x) 

X 



<CoQ 



n+l 



where 



2C2 - ci (ci + C2)^ ^ 3 _2 

Q ■= — — — < 1' aiid Co := — :^^^-^ — < -q 



2C2 + Ci 



2C]^C2 



(4.2) 



We remark that B is bandlimited, and that B-' (^, C) = for dist(^, C) > j^q- It follows that for ^, C 
with reFg < dist(C,^), we have pn(B){£,X) = 0- Thus we can use the spectral theorem for positive 
operators to write 



|B-i(e,C)| = |B-i(C,C) -Pn(B)(e,C)| < IIB-I -p„(B)||2 = max 

xeo-(B) 



1 



■Pn{x) 



<CoQ 



n+l 



By setting 



r := -logQ = -log 



we derive the following estimate on |B ^(i^,C)h 



|B-i(e,C)|<-cr2exp 



2C2 - Ci 
2C2 + Ci ' 



Tq 



We note that r depends only on ci and C2 and not at all on T. 

As we did for R, we estimate the L^o operator norm of B^^ to obtain 



|B-i Hoc < U'' E Kmik + l^e-^''/^ < er« 



where 



fc=o 



e = 3cr^KMT-'^ 



r'^e '"dr. 



(4.3) 



To complete the proof, we note that ||B ^R||oo < CI"^ x ^F '^S ^ < |, and consequently that 

||G-i||oo <2er'^. D 

Remark 4.2. In J^ \17}j it is shown that invertible matrices that decay exponentially away from 
the diagonal have inverses with the same property. However, it is not clear that their results are 
applicable to this situation. In particular, it is unclear how to control the ensuing constants in order 
to bound the norms of these inverses. 



5 Main Results 



In this section, we show that for a space X with a local, stable basis, the L2 projector is uniformly 
bounded in each Lp norm. 
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Theorem 5.1. Let H C M &e a quasiuniform set. Assume that h3. 1\) and As sumption \3.4\ hold 
true for the basis (i'^)^ge: of X. Then for all 1 < p < oo, the Lp operator norm of the L2 projector 
r= = V{V*V)~^V* : Li(M) -^ X is bounded by a constant depending only on M, p and the 
constants u, Ci, ci and c^ that appear in \3. i|) and Assumption\3.J\ 



Proof. We note that Proposition 13.61 provides upper bounds for the norms ||F||oo and ||V^*||oo) with 
which we can easily write down 

||F||oo = maxy^|f2,f(a;)| = q^'^^'^ maxy^ \voo,Ax)\ < q~'^^'^C2, 

X ^ ' X ^ ^ 

||y*||oo = max / |t;2,5(a;)|dx < g'*' ^ max / |^;i_^(x)|dx < g'^'^C2. 

It follows that the product ||l^||oo and ||I^*||oo is bounded by c^. In Proposition 14. H we have shown 
that ||(V^*V^)~"'^||oo is bounded. Together, these imply that 

||TT_|| IITI-II <^ IIT/II IIT/*II II ('T/*T/^~l II 

Thus by Riesz-Thorin interpolation, the Lp operator norm (1 < p < cxo) satisfies the same bound. 

D 

An immediate corollary is that for quasi-uniform sets H, the L2 projector Th gives rise to universal 
near-best approximations. 

Corollary 5.2. Under the conditions of the previous theorem, for f ^ Lp 

||/-Ts/||p<(l + 3^)distp(/,5(K,H)), 

where % < max||r='|L < 00. 
p 

A key motivation for studying stability of the Lagrange basis was the observation in [8] that the 
Lagrange functions associated with the Matern type kernels Km,M satisfied Assumption 13.41 In the 
previous section, it was observed that they satisfy the stability relation (|3.ip as well. It follows that 
Theorem 15.11 holds in the cases considered in p] . 

Corollary 5.3. Let Wl be a d- dimensional, compact, Riemannian manifold, let m > d/2 and let 
i^m,M be the kernels considered in 18]. If 3 cM. has mesh ratio p, then the L2 projector is bounded 
on Lp, 1 < p < 00 by a constant depending only on p, M and m, and, therefore delivers near-best 
approximation in each Lp. 



5.1 Example: L2 projection on §^ 

We now provide the example of the kernel K2 §2 that satisfies our two main assumptions, and 
which, by Theorem O gives rise to uniformly bounded L2 projectors for quasi-uniform centers. 
Furthermore, the Lp approximation orders for this kernel can be given exactly, which furnishes the 
operators T= with a precise error analysis. 
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These facts follow from an analysis carried out in [8] , where ^2 §2 was observed to satisfy Assumption 
3.41 a fact which was crucial to yielding bounded Lebesgue constants. A complementary fact, 
observed in [81 Section 5], is that K2 §2 belongs to the class of SBFs studied in [12j. A consequence 
of this, is that the kernel networks have universal convergence rates, a fact we consider below. 



Setup: Let Ld := y A^ — Agd be the pseudo-differential operator where Agd is the Laplace- 
Beltrami operator on S'^ and A^ = ^^. The Bessel potential Sobolev spaces have the norm given 
by 



m 



Y^i^ + XdVPif 



i=0 



LP{S<i) 



We note that when 7 G 2Z, L^ is a differential operator and Wp , the usual Sobolev space (imported 
from R via a partition of unity and local diffeomorphisms, or obtained using covariant derivatives) 
is embedded in H^ (this does not hold for other values of 7 when p = 1, 00). 

Finally for /? = 1, 2, . . ., let Gjs be the Green's function for L^, i.e., 

L^Gpix ■ y) = Sy{x). 

More generally let cpi^ := G^ + G^ * Tp where ip is an L^ zonal function and let 5= := span{(/)i^{x ■ 
S,): ^ € S} where S has uniformity measure h/q < cq. 

Convergence rates: The following was given in [12, Theorem 6.8] 

Theorem 5.4. Let 1 < p < 00, (3 > d/p' {1/p + 1/p' = 1). /// G H^, then with (pj3 and S~, as 
given above, 

distp(/,5H)<CA^II/llH- 

In particular, we have for /3 G 2Z that 

distp{f,SE)<Cp^h^\\f\\^,. 
by the embedding stated above. 

One can also treat functions of lower smoothness using a standard ET-functional approach. To this 
end, we assume /3 G 2Z, ensuring that Wp C H^. For < s < /3 and 1 < p < 00, we define the 
Besov space B^ ^ as the collection of functions in Lp for which the following expression 

\\f\\Bs^:=snpt~'K{f,t) 
is finite, where the iT-functional K{f, •) : (0, 00) — )• (0, cxo) is defined as 

i^(/,t):=inf{||/-<7||L,+t''||<7||^^^^: g^W^]. 

For a complete discussion of Besov spaces, we recommend [20 [ 1.11, and Chapter 7] and the ref- 
erences therein. An immediate consequence of the definition is that, for an arbitrary function in 
-Bpoo, and t > 0, there is gt G Wp simultaneously satisfying 



gt\\w0 < 2t'-"\\f\\Bs^^ and Wf-gth, < 2f\\fU 



p,oo 



In the setting of the previous theorem, one can approximate / by the kernel approximant to 5ft. 
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Corollary 5.5. Letl<p<oo, /3 > d/p', < s < /3, {1/p + 1/p' = 1). If f e B^,^, then with 
(j)p and S-E as given above, 

Example (The case m = 2, M = S^.) In [8i Section 5], the example K2S2 is considered. We note 
that 2 > 1.5 = d/2 + .5, so by Lemma 17.31 its shifts satisfy Assumption 13.51 with e = .5. 

By applying a simple symmetry argument, the kernel is shown to be a function of the the collatitude 

0E[O,7r], 

'i2,E.2ix,y) = (j)(0), where 6* := cos"^(x,y). 



I.e., it is zonal. Further analysis shows that it satisfies the fourth order differential equation 
^ (sin 96")" - -i- (cos 6 cot 96')' - -^ (sin 96')' + 6 = 0. 



sm U sm t) sm ( 

By considering a power series solution of this ODE, selecting the sole solution with appropriate 
singularity near 9 = 0, it is observed that (j) = G13 + Gp * ip with /3 = 4. Consequently, for 
quasiuniform centers the Lagrange functions associated with this kernel provide a stable basis, 
interpolation is stable in L^o, and L2 minimization is stable in Lp for every 1 < p < 00. Moreover, 
the abstract convergence rates from Theorem 15.41 and Corollary 15.51 combined with Corollary 15.21 
imply the following: 

Theorem 5.6. Let 1 < p < 00 and suppose H C S^ is quasiuniform, with mesh ratio p. If f £ H^, 
then 

\\f-TEf\\p<Ch''\\f\\H.. 

If0<s<4 and f e B^,^, then 

\\f-TEf\\p<Ch^f\\B^^^. 



6 Conclusion 



A major practical drawback in using the Lagrange basis is the computational overhead needed 
to construct it. A fundamental challenge will be to construct local kernel bases that satisfy (13. ip 
efficiently. Another, mainly theoretical, challenge is to find kernels satisfying the two basic assump- 
tions of this paper for which precise Lp error estimates are known, as was the case in the previous 
section for ^2 §2. A third issue is to find kernels having explicit, closed form representations (unlike 
the case K2 §2 , which is given as the solution of a fourth order ODE) . 

These last two points lead to the following conjectures 

1. Kernels on spheres: The first conjecture deals with approximation by restrictions to spheres 
of the surface splines. I.e., these are kernels used on R restricted to § . Earlier analytic 
properties of these have been studied in |10 1 ll3 t [2T]. The restricted surface splines are known to be 
conditionally positive definite. 
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For centers H C S , define kernel spaces 

generated by shifts of tlie restricted surface splines 

/(l-(x,y)r"'^/2log(l-(x,y)) deven 

'^"^^'^'^^■"\(l-(x,y)r"'^/2 dodd ' 

and spherical harmonics of degree \m — d/2\ or less. Such kernels are known to be of the form 
Gp + Gp * tp, with /3 = 2m |121 (3.6)] and thus have a precise, theoretical error analysis (similar to 
the case of K2 §2 of the previous section) . They also invert an elliptic differential operator of the 
form n^i(^ ~ O')' with ri, . . . ,rm a sequence of real numbers [TJ Lemma 3.4], as such, they are 
very similar to the kernels K^^d, which also invert elliptic differential operators of order 2m. 

Conjecture 1. If 'E is quasiuniform with mesh ratio p, the Lagrange functions Xi ^ S{ipra,'S) 
are a local basis (in the sense of Assumption \3.4\ ) and satisfy the Holder continuity assumption, 
Assumption \3.5[ 



The validity of this conjecture would have several immediate consequences: 

[A] The associated Lebesgue constants for interpolation are bounded by a constant depending 
only on p (and not on #S), 

[B] Interpolation provides precise L^o approximation rates. For functions f ^ G'^ , s < 2m, one 
has 

\\f-hf\\oo = 0{hn. 

[C] The Lagrange basis satisfies the stability condition (13. ip . 

[D] The Lp norms of the L2 projector Ts associated with these spaces are bounded with constants 
depending only on p. 

[E] The L2 projector provides precise Lp approximation rates. Le., approximation order s for 
functions having Lp smoothness s, for s < 2m (as determined by inclusion in Bessel potential 
i?2m ^^ Besov spaces Bp^): 

\\f-Tsf\\p = om 

A more ambitious problem would be to prove the previous for all kernels of the form 

(t>i3 = Gfs + Gp*ip, P> d, 
which includes the compactly supported SBFs, as well as most other SBFs of finite smoothness. 

Conjecture 2. //H is quasiuniform with mesh ratio p, the Lagrange functions Xi S S{(j)p,'S) are 
a local basis (in the sense of Assumption [3^ and satisfy Assumption \3.5[ 
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As a corollary, [A]-[E] would hold for kernels (pi^ as well. 

2. Surface splines on SO{3): The second conjecture deals with an analogous problem on SO{3). 
In this case, we consider the family of kernels 



km{x,y) := {sm{uj{y ^x)/2)) 



m-3/2 



where uj{z) is the rotational angle of z G SO{3) (the expression uj{y ^x) is the geodesic distance 
between rotations x and y). For centers H C SO{3), define kernel spaces 



jyl^mi ^) '■ — * 



Y,M^i-^^) + P\ P e nr„_3/2] and ^AgQ(0 = 0, VQ G Hp 



m-3/2] 



where H 1-^-3/2] is the space of Wigner L'-functions of degree \m — 3/2] or less. As in the previous 



case, they invert a differential operator of the form ]^"^]^(A— r^), [9] Lemma 3], and have theoretical. 



precise error estimates. 

Conjecture 3. /f H C SO{3) is quasiuniform with mesh ratio p, the Lagrange functions Xi ^ 
S{km,'^) are a local basis (in the sense of Assumption [XTJ ); and satisfy A ssumption 1 3. 51 

In this case, points [A]-[E] also would follow. 



7 Appendix 

In this section we relate the size of a Holder seminorm of a function to its norm in a Sobolev space. 
It is easy enough to state and prove in the Euclidean setting. 

Lemma 7.1. Let e > 0, m > d/2 + e, r > and suppose the finite set of points H C i?(0,r) is 
sufficiently dense so that h := max2.g^(o,.\ dist(x,H) < hor (with Hq := l/{32m'^)). Then there is a 
constant so that for every f S W^(^B{0, r)) and y G B{0, r) 

|/(y)-/(0)|<C|yrr-"^/2-e| ^ f \Dy(x)fdx] . 
Proof. Suppose r = 1. By the Sobolev embedding theorem, there is a constant so that 

\u{Y)-um ^^/ ^ /■ ,^^ ^^. 



\Y[ 



<C I y /" \D''u{X)fdX 



|a|<-m ■ 

for all u € W^{B{0, 1)). By rescaling this with rY = y and f{y) = u(Y), we have 

1/2 



l/(y)-/(o)l 

\y/r\^ 



<C7 I V r^l^l-"' f \D''f{x)fdx 

\ , T^._ JBiO.r) 



1^ JB(0,r) 

\a\<m ^ ' ' 
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for all / S W^{B{0,r)). We can now apply [U Lemma 3.7] (a standard Bramble-Hilbert type 
argument) to obtain 

C I V r^l")-'* / \D''f{x)\^dx] < Cr""-'^/2 [ y- f \D''f{x)\'^dx] 

\ I T::^ JB(0,r) / III JB(0,r) } 

\\a\<m \ ' / I \^a\=m y ' ' J 

which proves the lemma. D 
Before proceeding, we need some necessary geometric background, which can be found in [S]. 



Exponential map We define the exponential map Exp„: M -^ M, our basic diffeomorphism, by 
letting Expp(O) = p and Expp(stp) = 7p(s), where jp{s) is the unique geodesic that passes through 
p for s = and has a tangent vector 7p(0) = tp of length 1. 



Metric Equivalence The family of exponential maps are uniformly isomorphic; i.e., there are 
constants < Fi < r2 < oo so that for every po gM. and every x, y G -6(0, tm) 

ri|x-y| < dist(Exppg(3;),Expp^(2/)) <r2|x-y|. (7.1) 



Sobolev spaces Let il C M be a measurable subset. We define the Sobolev space W2^{Q) to be 
the space of all functions / : M ^^ M so that the real valued function p i-> [ V'^/|c,,p = (V™/, '^^f)g,p 
defined for p G Vt \s m. L2. The associated norms are: 



m „ \ 1/2 



X; /ivVlLdMp) ■ (7.2) 

7 — n .^ ^ / 



\A:=0" 

For the precise meaning of (V™"/, V'"/)g,p , we refer the reader to jHJ Section 2], although we note 
that in R'^, its exact meaning is 



E («) ^"/(P)^"5(P), 



\a\=k 

i.e, the inner product of the fcth order derivative tensor at p, in which case (|7.2|) is the usual 
Euclidean Sobolev norm. 



Metric equivalence for Sobolev spaces The exponential map also gives rise to the following 
metric equivalence, proved in [8, Lemma], between Sobolev spaces on regions in the manifold and 
Sobolev spaces on regions in W^. For ?tt, G N and < r < tm, there are constants < c < C so that 
for any measurable Q C Br, for all j € N, j < m, and for any po € M, 

c\\u o Expp^ ||^^,(^^ < ||n||^.(E,p^^(^)) < C\\u o Expp^ H^^.^^^ (7.3) 

We are now able to state the analogous zeros lemma for manifolds: 
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Lemma 7.2 (Zeros Lemma). Let e > 0, let m be a positive integer, greater than d/2 + e, and let r 

be a positive real number less than tm, the injectivity radius o/M. Suppose that S C B{x,r) C M is 
a discrete set with h < Tihor. If u G W^{B{x,r)) satisfies u\^ = 0, then for every x,y G B{x,r), 

\u{x) -u{z)\ < Cr™"'"'^/^dist(2;,z)'||n||vi/m(B(x,r)), 
where C is a constant independent of u, x, z, h and r, and where Pi depends only on m and M. 



Proof. This follows from (j7.ip . (j7.3p and Lemma |7.H (which we can apply since h < ^j^r implies 
that the J^xp^(3) has fill distance less than h^r) in precisely the same way as |8, Lemma 3.8]. 
Letting z = Exp^(y) 

\uix)-u{z)\ = \u{0)-u{y)\ < r'^-'^/^-'/^\y\'\\u\\w^^B{o,r)) < Cr'^-'^/^-'[d^st{x, z)Y\\u\\w.Y^^B{x,r))- 



We note that \y\ = \y — 0\ = dist(Exp^.(y),Exp^(0)) = dist(z,x). 



D 



We now show that the Lagrange functions associated with the kernels Km,M satisfy the Holder 
continuity property, Assumption 13.51 

Lemma 7.3. Let < e < 1. For m > d/2 + e and for H C M having mesh ratio p, there is a 
constant C depending only on m, M, e and p for which 



\xdx)-xM\<C 



dist(j;, y) 



Proof of Lemma \7.3[ By compactness and continuity of the functions s, it suffices to consider cen- 
ters H and pairs x,y that satisfy dist(x,y) < h < Tho, a quantity appearing in the zeros lemma, 
Lemma 17.21 (hn is defined in Lemma |7. II and Pi is the lower bound in ()7.ip ). 

By [SI Corollary 4.4] there is a constant zv > such that, then for < T < rM we have the estimate 

where q is the minimal separation distance between points of H. We proceed as in the proof of [H 
Proposition 4.5] by considering three cases. 

Case 1. For dist(^, x) < ^-y~! we can apply Lemma[72]with r = 2j^^ (i.e., on the ball B{x, 2^^^), 
where the zeros of x^ have mesh norm bounded by 2h) to achieve 



xd^) -xdv) 



[dist(x,y)]^ 



<C 2 



/iqPi 



m—e—d/2 



\X^\\m,B(x,2j-^) - ^ 



m—e—d/2 



dist(5,a;) n 



Thus, for X G B{i, ^^) we have \xi{x) - Xi{y)\ < C[dist{x,y)Yq 'exp(-z^ — j^ — ;. 
Case 2. A similar argument holds for ^y- <T = dist(^,2;) < rM, where we obtain 

\x,{-)-xdy)\ < c^^™-^-'/^[dist(x,,)]lxdl^^.(M\B(,,T-4.)) 

/dist(x,^)~ 



< C[dist(x,y)]'^g ^ exp 



I- 
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In the second inequality we apply [U Corollary 4.4] after applying Lemma 17.21 and making the 
observation that B{x, ^) C M \ B{C, T - ^). 

Case 3. For T > rM, we make use of the fact that M \ B{S,,icm) D M \ B{S^,T), so we simply use 
the estimate for T = rm to obtain 



"' h 



\xd^) - Xdy)\ ^ C[d\si{x,yyfq 'exp 
This completes the proof of the lemma. D 
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